Optical Lattice Hamiltonians for Relativistic Quantum Electrodynamics 



O 

> 
O 



cr: 



I 

o 



O 



X 



Eliot KapiiE and Erich J. MucUeiQ 

Laboratory of Atomic and Solid State Physics, Cornell University 

We show how interpenetrating optical lattices containing Bose-Fermi mixtures can be constructed 
to emulate the thermodynamics of quantum electrodynamics (QED). We present models of neutral 
atoms on lattices in 1+1, 2+1 and 3+1 dimensions whose low energy effective action reduces to that 
of photons coupled to Dirac fermions of the corresponding dimensionality. We give special attention 
to 2+1 dimensional electrodynamics (QED3) and discuss how two of its most interesting features, 
chiral symmetry breaking and Chern-Simons physics, could be observed experimentally. 
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I. INTRODUCTION 



Planar quantum electrodynamics (QED3) has been of 
great theoretical interest for many years. The restriction 
to lower dimensionality profoundly changes the structure 
of the electromagnetic field and its symmetries [l| . In two 
spatial dimensions (2+ld) the electromagnetic potential 
between point charges that scales as logr instead of r~^, 
leading one to expect charge confinement in the appro- 
priate limit. QED3 with massive fermions is particularly 
rich. One of the two possible Lorentz invariant mass 
terms breaks parity and time reversal symmetry, leading 
to a topological Chern-Simons term in the photon ac- 
tion. This term produces a gap in the photon spectrum 
and makes the effective interaction between the fermions 
short ranged. The other Lorentz invariant fermion mass 
term breaks "chiral" symmetry and is spontaneously gen- 
erated at low temperatues. Both confinement and chiral 
symmetry breaking occur in quantum chromodynamics; 
since QED3 is a much simpler theory it has been studied 
as a model system for these effects 043- Further, QED3 
has been proposed as an effective theory for spin-liquids 
d. 111 and the pseudogap phase of high-temperature su- 
perconductors [lOl4l3l | . A direct experimental realization 
of QED3 could therefore be relevant to many areas of 
both condensed matter and high energy physics. 

Recent advances in cold atom physics and optical lat- 
tices [l3| have made the simulation of QED3 a real ex- 
perimental possibility. Dirac fermions app ear at half 
filling in a hexagonal optical lattice |15l4l9l |. mimicking 
the relativistic band structure of graphene [23|. They 
also appear in a square lattice with appropriate complex 
hopping matrix elements [2]| . This has led to a num- 
ber of proposals to emulate relativistic physics with cold 
atoms |2l| - |2a |. Here we extend these previous proposals 
by showing how to emulate QED3 with cold atoms. By 
combining mixtures of bosons and fermions on interpene- 
trating optical lattices with appropriately tuned densities 
and interactions, both relativistic fermions and the full 
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3-component dynamical gauge field can be constructed, 
allowing QED3 to be probed directly. The linear Bo- 
goliubov excitations of three species of bosonic atoms 
correspond to the three components of the vector field 
Af^. One of the boson species occupies the same lattice 
as the fermions, while the other two sit on alternating 
bonds between the fermion sites. These bosons modu- 
late the fermion hopping amplitudes, leading to an inter- 
action that at low energy has the form of a gauge cou- 
pling eijj (7^j4p) ijj. Being a lattice model, higher energy 
terms in this theory will naturally deviate from the emer- 
gent Lorentz and gauge invariances, but at low energies, 
these terms are heavily suppressed, and after applying 
the Wilsonian prescription of discarding irrelevant oper- 
ators, the resulting Hamiltonian is that of QED3. 

It is important to note that in the real time formalism, 
the time {Aq) component of the gauge field is canonically 
quantized by operators with negative norm, with causal- 
ity and Hermiticity preserved in real states by the Ward 
identity [l^. This is a general feature of vector bosons, 
and cannot be duplicated with cold atoms. However, in 
the finite temperature (imaginary time) formalism, all 
the components of the A field have positive norm. Thus, 
we focus on producing a quantum simulator which repro- 
duces the thermodynamics of QED3. Many of the most 
interesting properties of QED3 (such as mass generation 
and symmetry breaking) can be probed through equilib- 
rium thermodynamics. Further, some real-time correla- 
tions can be accessed through the fluctuation-dissipation 
theorem. 

Our construction is very different from traditional ex- 
amples of gauge fields in condensed matter systems 0- 
[ll,l23|- Typically, these gauge theories arise from using 
descriptions with redundant degrees of freedom (such as 
slave bosons in the large U Fermi Hubbard model). The 
complicated underlying physics of these emergent gauge 
theories makes it difficult to use them as simulators of 
gauge physics itself. The approach that we take is differ- 
ent; instead of formulating a redundant description of a 
strongly interacting system, we take a weakly interacting 
system with a large number of degrees of freedom and 
tune the couplings and particle densities so that the low 
energy physics is identical to thermal QED. 

To construct our theory, we necessarily choose a par- 
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ticular gauge. The Hamiltonian of our theory is identical 
to that of QED in the Feynman gauge. Our theory is 
therefore not a true gauge theory in the formal sense of 
being described by an equivalence class of distinct actions 
and wavefunctions; rather, it is a theory of a relativis- 
tic fermion interacting with three neutral, independent 
massless bosons through a coupling of the form ip^f^A^ip. 
All physical observables of QED are gauge independent, 
and despite its fixed gauge, our simulator should correctly 
reproduce all such observables. 

The remainder of the paper is organized in the follow- 
ing manner. In section we will describe the canoni- 
cal quantization of QED3 in the Feynman gauge, taking 
care to denote the various terms in the Hamiltonian that 
must be replicated in our cold atom system. In section 
mil we will then outline our model and show that its low- 
energy limit is identical to the Hamiltonian formulation 
of QED3. At the end of that section we will describe some 
of the more dramatic properties of QED3 and how they 
would manifest themselves in cold atom experiments. Fi- 
nally, in sections IIVI through IVII we will offer concluding 
remarks and discuss extensions of our theory to 1-1-1 and 
3+1 dimensions. 



II. QED AND GAUGE FREEDOM IN 2+1 
DIMENSIONS 



Our goal is to simulate the Euclidean action for 2+1 
dimensional quantum electrodynamics, 

Le = iijj'-fitdfj.ip + ei' (7m^m) - ^F^j,yF^^. (1) 

We adopt the conventions that spacetime indicies which 
run from to 2 arc represented by Greek letters, and 
objects in boldface will represent the two spatial com- 
ponents of the Lorentz 3-vectors. Our imaginary time 
metric and 7 matricies are 



70 = «cr2,7i = ^(^x,l2 



(2) 



and ■0 is a two-component fermionic spinor. The Uxyz are 
Pauli matricies, and F^i, = d^j^A^ — duA^ is the electro- 
magnetic field strength tensor. Summation is implied for 
all repeated indicies, and we have abandoned the stan- 
dard convention of raising and lowering indicies to em- 
phasize that our imaginary time metric is trivial. In sec- 
tion HI we will generalize this action to massive fermions. 
We work in the imaginary time formalism [2^, [s^ to en- 
sure that the photon propagator D^^ is positive definite 
and therefore all photon states have positive norm and 
there are no ghostlike degrees of freedom. We note that 
in a non-Abelian theory, negative norm degrees of free- 
dom will be found even in the imaginary time formal- 
ism. Thus, our approach is not appropriate for simulat- 
ing QCD. 



We will simulate this Hamiltonian by designing an op- 
tical lattice system that exhibits identical degrees of free- 
dom and Feynman rules up to lowest order in powers of 
momenta. To arrive at a simple enough form of the ac- 
tion for quantum simulation, we will employ the Faddeev- 
Popov gauge fixing procedure . This scheme enforces 
the gauge constraint 



w (x) 



(3) 



and then path integrates over w to remove the gauge de- 
grees of freedom from the Lagrangian. So long as we re- 
strict ourselves to calculating only gauge-invariant quan- 
tities, the integration over w can be performed trivially. 
This alters the Euclidean action by adding a term of the 
form 
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Le ^ Le - — (5^^^.) : 



(4) 



where ^ is an arbitrary parameter. If we then choose 
^ = 1 (the Feynman gauge), the off-diagonal terms in 
the photon propagator cancel out and we are left with 
the action 



Lf 
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iilf,df,ip + ei (7^^^) ip + -df^A^d^A 



(5) 



Within this gauge, our system is described by the Hamil- 
tonian 



H = VFip{~<^ + ei{azAQ + ■ A)i 



(6) 
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-d^A^d^A^ = X! 5Z '^kaj^k 
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The Ur operators annihilate a photon of polarization r, 
and uj\!, — Vs\k\] this is simply the canonical quantization 
of the free photon Hamiltonian. Note that, unlike in 
Minkowski space, [a^k, ajp] ~ 5rs5^ (k — p) describes 
only positive norm states, since the negative time-time 
component of the photon propagator is made positive by 
the rotation to imaginary time. The gauge field A is 
quantized in terms of these bosons by 
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^ — ik-: 



Vk^ 



(7) 



Here, the /rk/i are real, 3-element polarization vectors 
satisfying /^-k • /sk = <5rs ■ The Ward identity [sij implies 
that real photons are transverse, ie k • f^k = for all 
physical states; only transverse photons will be produced 
or destroyed in observable processes. The momentum- 
space propagators for the fermion and gauge field are 

Dp = ,„ 7^ , . , (8) 



m + le 
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We will now show how cold atoms can be used to directly 
implement the Hamiltonian 
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FIG. 1: (Color online) Square lattice for QED3. The fermions 
are confined to the lattice sites K and S denoted by white and 
black circles. The bosons sit on the red triangles and blue 
squares denoted by a\ and Q2. Overall phases for fermion 
hopping are labeled in the figure. 

III. 2+1 DIMENSIONAL OPTICAL LATTICE 
QED 

We give two alternative geometries for the lattice con- 
fining our fermions. We require tliat tfie system exhibits 
"Dirac points" instead of a Fermi surface; all zero- or 
low-energy fermionic excitations must occur near some 
finite set of points K in momentum space, and have 
the dispersion relation ep = ±^yvjpp^ + m^ for small 
p measured with respect to K. The most well-known 
system of this type is the hexagonal lattice [l5l - [l9j , re- 
alized in condensed matter in graphene [20j and carbon 
nanotubes. In cold atoms, the ability to directly manipu- 
late the phases of lattice hopping amplitudes enables the 
creation of Dirac points in a square lattice [2l|, , some- 
thing that hasn't been realized in materials. We describe 
both approaches. 

A. Square Lattice Implementation 

Our square lattice approach is based on the geometry 
proposed by Liu et al [2l| to study the anomalous quan- 
tum Hall effect. Consider noninteracting fermions in an 
square optical lattice with nearest neighbor hopping t 
and lattice spacing /. We add to the hopping amplitude 
a set of phases such that 

H = (9) 

tjk = t {j/j = j/fc} (horizontal hop), 

= iti-lf'+y^ {yj- 7^ yfc} (vertical hop). 



The operators ijjj and '0j destroy or create a fermion at 
site j. This structure is mathematically identical to a 
constant magnetic field with flux density tt (or half of a 
fundamental magnetic flux quantum ^o) per plaquette. 
This effective magnetic field could be realized through a 
spati ally dep endent Raman coupling to an external laser 
field [Hif . 

To find the band structure, we break the lattice into 
two sublattices R and S, depending on whether Xj + yj is 
even or odd (figure [T|) . The momentum space Hamilto- 
nian can be written as the dot product of a pair of two- 
component spinors through a matrix. This matrix M has 
no diagonal elements because nearest-neighbor hopping 
always transfers a fermion between sublattices. 



Hd = -2t^{^k'V'L} -M.H'Bk^^Sk}, (10) 

k 

_ / COS k^l + i COS kyl \ 

\ cos kxl — i cos kyl J ' 



This Hamiltonian produces a Dirac cone dispersion about 
the four K points K±± {±7r/2Z, ±7r/2/}. Only two of 

these points (K_|--|_ and K_| ) arc distinct; the others are 

related to them by symmetry. We predominantly focus 
on K++ = 7r/2Z{l,l} = K and let k = K + p; the 
structure at the other K point is the same up to the 
signs of the Dirac 7 matricics. With the definitions ([2|), 
small p and tp-p — {V'flp, "^Spli our Hamiltonian becomes 



= 2«J]*p(T-p)?'>p, 

P 



We now add a gauge field to the theory by introducing 
three species of lattice bosons. One of the three boson 
species (ao) lives on the same lattice as the fermions and 
the other two (ai and 02) occupy interpenetrating square 
lattices as shown in figure [TJ If the interaction between 
the bosons and fermions is repulsive, the presence of oi 
bosons between a pair of neighboring fermion lattice sites 
acts as an effective potential barrier that reduces the hop- 
ping amplitude between the two sites. Similarly, the pres- 
ence of the fermions will modulate the boson hopping and 
lead to a nearest-neighbor repulsion term between bosons 
and fermions. The total Hamiltonian incorporating all of 
these effects is 
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{x e Lf, yb=i,2 e ib=i,2} 



(12) 



x,s 6;k x;6 X 

+ 3''V'lj,/Sx^S/i?,x±s"by.=x±f +Xl5b'^6y.V'lj,/SxV'i?ySx+ ^1 3&4y.±2s«byb V'l?ySx=y,±f V'i?/Sx=y.±l 

b=l,2;y 



6=l,2;x,s 



b;x 



The sum over b is over boson flavors (either 0,1,2 or 1,2 
as specified) , tlie vector s joins nearest neighbors and is 
equal to {±lx,0) , (0, ±Zy), where, as previously defined, 
I is the lattice spacing, i/ is the fermion lattice, and 
Lf, is the lattice associated with species b. The operator 
TT-by measures the density of ai, particles at y. From left 
to right, the terms on the first line are the fermion hop- 
ping term, the boson kinetic terms ef,k (which include the 
boson chemical potentials), the weak repulsive contact 
interaction Vf, between bosons, and the local repulsion 
between fermions and the ao bosons. The terms on the 
second line are the boson mediation of fermion hopping, 
nearest neighbor Bose-Fermi repulsion, and fermion me- 
diation of boson hopping. We assume that there are no 
direct interactions between bosons of different species as 
they are spatially separated. The gb are the effective in- 
teractions between the bosons and fermions, which can 
be tuned by manipulating the optical potential depth or 
through Feshbach resonances and are calculated from the 
overlap of the neighboring Wannier functions (36j . For 
simplicity, we shall assume that ffb = ff^, a- situation that 
we expect can be reached through appropriate fine tun- 
ing. In the tight-binding limit, next nearest neighbor 
hopping and any additional interactions are vanishingly 
small. The last two terms in Eq. can lead to un- 

wanted terms in our low energy action. As calculated be- 
low, their contributions cancel each other if the 6 = 1,2 
bosons are condensed about K = {it/21, tt/21}. This can 
be engineered, for example, by using the techniques of 
Ref. |33| to introduce phases on the hopping matrix ele- 
ments. Under these conditions, 



eb=l,2k 



eok 



+24=1.2 1 cos — 



cos 



(13) 



—2tQ (cos ^ 



' cos ■ 



where the tf, are positive numbers parametrizing the bo- 
son hopping. We shall now show that, with appropriate 
choices for the various parameters, the low energy exci- 
tations of eq. ([T^ are exactly those of QED3. 

We first assume that the temperature is low enough 
that all of the boson fields have condensed; if this is 
the case, we can expand the Bose operators around the 
condensate expectation values, dropping terms that are 
quadratic or higher in boson fiuctuations. For example. 



if 6 = 1,2 



^bx = 7^b + \/nb(-l) 

= nb + V^(-l)"+-^a 



O I 1(14) 



bx, 



where at,x — a-bx — ^/nt- To treat the bosons' self- 
interactions and vacuum expectation value, we perform a 
Bogoliubov transformation [s^]- For a given species this 
process gives 



Hb = Ebual^abk, 

k 

Ebk = \/ebk (ebk + 2nbT4), 



flbk 



Mbk^bk - Wbka^_k 



Here Ubk and Vbk are real constants satisfying u^ — v^ = 1, 



"bk = 



Vbk 



Cbk + nbVb , 1 



bk 



2' 



(15) 



/ £bk + nbVb _ 1 



Ekt 



Furthermore, at low momenta these coherence factors 
simplify. 



lim Wk 

k^O 



lim Vk = \/nV/Ei^, 

k^O 



allowing us to use the low-energy approximation 

Gbk 



2\/n,£ 



abk + oib-k 



(16) 



(17) 



The "speed of light" for these bosons is Vb ~ y/nbVb/mb. 
We shall tunc the tb so that the Fermi velocity vp 
matches each of the three Vb- 

Now consider the first term on the second line of Eq. 
(|12p . Expanding the fields in momentum space and sum- 
ming over k and b yeilds 



Y 5b'0fl/SxV'5/flx±l"bx±li,/2 



(18) 



b=l,2:: 



E 

b=l,2,k,p,q 

■ Y 9bnb (■0k7bV': 

b=l,2;k 



ribgb (Gbk'0p7bV'q ) (^k+p-q 



k J + O ( -,/k, /p, ?q 



5 



Here, the momenta in G and ij: are expanded about 
{11/21,71/21} and K, respectively. Similarly, the local re- 
pulsion term adds two terms of the form 

E (V'kTo^k) (19) 

X k 

+,gO\/nO ^ (^GokV'p70'/'q) (^k+p-q 
k,p,q 

+0 f-,/k,?p,/q) . 



We can account for the remaining two terms on the sec- 
ond line of (|12p through a similar expansion. Assum- 
ing that Qb = g'b, we find that with our choice of boson 
dispersion (with a minumum at {7r/2/, 7r/2/}) these two 
terms precisely cancel each other to lowest order in k. 
No such cancellation occurs if the bosons are condensed 
about k = 0. The resulting low energy theory is 



H = vf'^ V'k7 • k^Ak + X! X! '^•'^"fek^bk + X! X! (V^Gbk + m) '(/'pTbV'k+p + O ( - , ^k, Ip, /q 



h=0 k 



6=0 k,p 



(20) 



Comparing Gbk in ()17p to the canonical quantization of 
the gauge field A in (O leads us to write 



Gh 



2V2^Ab 



(21) 



Provided that = ■\/ ni,Vb/mi, ~ vp and 2nb\/2Vbgb = e 
for all three boson species, we finally see that ((20)) reduces 
to ^ with the addition of a constant shift of the gauge 
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(22) 



This constant shift has no effect on the overall physics 
of the theory and simply shifts the location of the Dirac 
points in k space. 

For K^^ fermions and RW fermions with a lattice 
depth of 4£'fl (where Eji is the recoil energy) for each lat- 
tice, an average boson density of unity, and a Bose-Fermi 
scattering length of — 177ao [l^, the coupling constant 

e (in units of square root temperature) is approximately 

1/2 

2(nK)''". This value can be tuned considerably through 
Feshbach resonances, optical lattice tuning, and adjust- 
ing the boson filling. 



B. Hexagonal lattice 

Our theory can also be realized on a hexagonal optical 
lattice |15M20l|. Assume that the spacing between verti- 
cies in the hexagons is /. The basis vectors b connecting 
a pair of nearest neighbor sites are 



I 



2V3' 



After dividing the lattice into two triangular sublattices 
K and 5* (see figure [2]) we arrive at the fermion hopping 
Hamiltonian 



k j = l 



'''■''V'kV'5k + /i.C. 



(24) 




FIG. 2: (Color online) Hexagonal optical lattice for QED3. 
The fermions and ao bosons sit on the verticies of the hexag- 
onal lattice (in black). The two other boson species a\ and 
are on the interpenetrating orthorhombic lattices; ai bosons 
are confined to the blue squares, 02 bosons are confined to 
the red triangle), and the green circles indicate sites shared 
by the two lattices. The hopping parameters of these lattices 
can be tuned so that the low-energy dispersion is simply 
k^/2m* for each of the three species. 



The coefficients e*'^ ''j vanish at six points K, of 

which only two are unique ~ the other four are related 
by reciprocal lattice vectors. We shall choose K± = 



± 



(0, 1). For small k defined relative to one of these 



two K points, the Hamiltonian is the familiar Dirac cone: 



(25) 
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As in the previous section, i/i = {ipR,ips} is a two- 
component spinor and ip = '0^70 • Wc again construct 
the gauge field by adding three BEC fields, as shown in 
figure (HI). The boson qq sits on the same lattice as the 
fermions, but oi and 02 are instead on the links between 
the fermion lattice sites. Specifically, ai is chosen to oc- 
cupy sites separated by 61 and S2 and 02 occupies sites 
separated by Si and 63. This structure breaks the Cg 



The couplings g are determined from the fermion-boson 
scattering amplitudes and from the geometric position 
of the bosons relative to the fermion lattice sites [s^ . 
Repeating the steps in equations (|14j|20p . and setting 

552 = = ^ (1 + ^) 5- (27) 

yeilds the QED Hamiltonian ([6|) with an effective charge 
e = 2\/3V ng and a speed of light vp = tl\/3/2. 



C. Observable Consequences 

1. Physical and Virtual Photons 

An important property of QED3 is the fact that there 
is only a single physical polarization for the photon. This 
means that in any experiment neither time-like nor lon- 
gitudinal photons can be created or destroyed. Conse- 
quently, in our emulation of QED3, no perturbation of 
the fermions will produce these "unphysical" photons. 
This is a dramatic observable. 



2. Spontaneous Symmetry Breaking 

At sufficiently low temperatures, thermal QED3 is ex- 
pected to undergo a phase transition to a state with bro- 
ken chiral symmetry and a finite fermion mass [11-0, Isol - 
|4]| . Chiral symmetry relates the two Dirac points of our 
lattice model, and as we will show below, chiral sym- 
metry breaking on our square lattice corresponds to a 
commensurate {tt, tt} charge density wave. 

To define the Chiral symmetry, wc combine the 
fermions at the two Dirac points into a four-component 



symmetry of the lattice, but if the coefficients are cho- 
sen properly the low energy theory will be rotationally 
invariant. As in the square lattice case, we require that 
the hopping matrix elements for bi and 62 to be opposite 
in sign to those of 60 • This will again allow the unwanted 
nearest neighbor repulsion and fermion mediation of bo- 
son hopping to cancel each other. Ignoring these terms, 
we're left with: 



(26) 

I 

Dirac spinor. Defining k = K± + p, we write 

H = z^(r.p),V (28) 

Ipp = {V'i?^K++p, V'SK++p,V'flK_-(-p,'0SK_+p} , 




The photon coupling becomes etpTfj,A^ip, and we have 
simply obtained a larger representation of the Dirac al- 
gebra. This theory is invarient under transformations of 
the form 

•0 — >■ exp (ieiCi) exp (ie2C2) ■0) (29) 

for any real constants ei and 62 and 4x4 chiral matricies 
Ci and C2 defined by 

These matricies anticommute with each other and with 
all three F matricies. They are the generators of a sym- 
metry transformation that simultaneously flips between 
the two sublattices and the two K± points. 

The low temperature phase of QED3 breaks this sym- 
metry, with (^f. tpk'4'k) 7^ 0. In our real space lattice this 
order parameter is a charge density wave of the form 

(Y^^k^i^) = E ^Lv^sx - E ^Lv^flx. (31) 

fe xGLs xGLr 

This order can be observed using light scattering [i^l 
or through single-site imaging of the fermions. Unfor- 
tunately the transition temperature for this spontaneous 
symmetry breaking is Tc ~ 10~^e^ [SE^. Using our pre- 
vious estimate e ^ 3(nK)^''^, we find that ^ lOOpK is 
beyond the reach of current experiments. 



J 
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FIG. 3: Hopping terms that lead to a chiral mass He- The 
phase factors of (±1) are labeled in the figure; the phase of a 
hopping term changes sign when rotated by 7r/4 and between 
sublattices R and S for the same direction. These phases 
correspond to a virtual magnetic flux of tt per plaquette (half a 
flux quantum). Since this virtual flux was already introduced 
to create the Dirac fermions, a chiral mass could be introduced 
simply through allowing next nearest neighbor hopping. 




FIG. 4: The signs of the hopping amplitudes in Eq. (|39p 
for a given sublattice in the hexagonal lattice. All signs are 
reversed for hopping within the other sublattice. Properly 
implemented, this term yeilds a chiral mass and will create 
gaps in the spectra of the 3 BEG fields. 



3. Mass terms 



In 2+1 dimensions, there are two linearly independent, 
Lorentz-invariant terms which lead to a massive fermion 



dispersion relation uj = iJ^nk^ + ttiq + r 

Hm = V-'kV'k, 
k 

He = tTOcX]^kClC2V'k- 



(32) 



The first, H.,n, breaks chiral symmetry, and can be ex- 
perimentally engineered by introducing a superlattice po- 
tential. In terms of the original fermionic sublattices R 
and S, we can write Hrr, as 



.xeL£ 



xeLn 



(33) 



This term leaves the photon modes qualitatively un- 
changed, and as explained in Sec. lIII (T2] generates a den- 
sity wave. 

The chiral mass term He is more interesting, as it has 
dramatic consequences for the photons. As explained 
below, it can be engineered by introducing longer range 
fermion hopping terms with appropriately chosen phases. 

In the presence of chiral-massive fermions, the resum- 
mation of 1-loop corrections to the photon propagator 
introduces a Chern-Simons term [l[ in the action of the 
form 



A 

Lcs = -^^isFa^eaisj, 



(34) 



with A ^ nic- This term is not found in 3+1 dimensional 
theories. Within the Feynman gauge it yields a modified 
photon propagator: 



-1 



'til' 



(35) 



with a photon mass gap |A|. The action Lcs is only 
gauge invariant up to a boundary term, reducing the 
symmetry of the system to a subset of all possible gauge 
transformations. The properties of QED3 with a chiral 
mass term are therefore sensitive to the topology of the 
2+1 dimensional space. 

In terms of the two sublattice fermion operators, the 
chiral mass term He is 



He = (V'ljK++pV'flK++p-V'sK++pV'SK++p 

p 

- V'kK- -f-p V'flK_ +p + +p V'SK_ +p) • (36) 



He is invariant under a transformation which simulta- 
neously switches the two sublattices and rotates between 
the two K points, but is not invariant under either trans- 
formation individually. The simplest perturbation which 
obeys this symmetry is is a next-nearcst-neighbor hop- 
ping term with direction-dependent phases, shown in fig- 
ure [3] The Hamiltonian for this term is: 
7nc v-^ 



H. 



NNN 



4 ^ 



Sx+Sj 



(37) 
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where the four Sj vectors join a site and its next near- 
est neighbors. The phases of these hopping terms change 
sign under a 7r/4 rotation and between the two sublattices 
R and S. Such hopping phases are natural in our model, 
Eq. p2|) , where the fcrmions experience an effective mag- 
netic flux of TT per plaquette: the phases correspond to 
a flux of 7r/2 through each triangular half-plaquettc. To 
lowest order in p, Hnnn reduces to He- Tuning the lat- 
tice to allow nearest neighbor hopping would therefore 
introduce a chiral mass gap into the model. 

Engineering this chiral mass in our optical lattice 
would allow us to study Maxwell-Chern-Simons electro- 
dynamics [ll-l3l. l4i - |46t . As seen in Eq. there will be 
a gap in the photon spectrum. Since the photon in our 
model can be identified with excitations about a BEC, 
this photon mass implies the existence of a neutral Bose 
condensate with a gapped spectrum. This unusual spec- 
tral feature could be detected through RF spectroscopy 
[4?! or through its thermodynamic consequences. 

Further, given that p in Eq. p6|) is small compared 
to K±, one can write, in terms of the d-density wave 
operator &: 

He - 2m,^ V'li/sk-H{^/,,^/,}^i?/SkSin26lk = 2m^(d 

R.S 

(38) 

where tan(0k) ~ kyjk^- Several authors, including 
Chakravarty et al [48|, have proposed that such d- 
density wave order plays a role in the pseudogap phase 
of high temperature superconductors. 

Connection to the Haldane Model 

On the hexagonal lattice the hopping term, analogous 
to Eq. ([57]) . leading to a chiral mass is (Fig. S]): 

E (-l)'(V'k^«r.+r,-V'L^5r.+r,).(39) 

The sum on k runs over all lattice sites in a given sub- 
lattice and Tj are the six vectors which join a site and 
its next nearest neighbors. This construction is analo- 
gous to the Haldane model [i^ for realizing conductance 
quantization in the absence of an external magnetic field 
in a hexagonal lattice. The contribution to the low en- 
ergy theory from this mass is the same as in the square 
lattice. In an experiment the appropriate phases could 
be engineered through the use of light assisted hopping, 
allowing the hopping elements to pick up local phases 
from an external laser beam 1351. 



4- Robustness 

An important issue with our proposal is that of sta- 
bility: how robust is our construction against symme- 
try breaking perturbations? Inevitably, any realiza- 
tion of QED through this or a similar scheme must 



contend with anisotropics in the couplings, unequal or 
anisotropic Fermi and Bose velocities, and residual in- 
teractions which violate our emergent gauge invariance. 
Further, we have ignored the effect of the harmonic trap- 
ping potential required to confine the atoms; the long- 
ranged forces of QED could be disrupted if the trap were 
not shallow. 

These are difficult issues, and we are unable at this 
point to quantify all of them. Each of these effects would 
have to be treated in detail, and their importance will 
inevitably hinge on technical details in the experiment. 
Some particular perturbations have been studied. For 
example, we have already seen that at low temperatures, 
QED3 is unstable towards spontaneous chiral symmetry 
breaking. Near the transition temperature, small pertur- 
bations which break chiral symmetry can be quite im- 
portant. Other perturbations are less destructive: for 
example, Vafek et al. [l^l , showed that anisotropy in the 
Fermi velocity is an RG-irrelevant perturbation at low 
energies. 

We do have some additional knowledge of the ro- 
bustness of this model agains Lorentz invariant per- 
turbations. Through a simple power counting, we 
argue that any generic Lorentz invariant perturba- 
tion will have a vanishingly small contribution to 
the low energy physics of our system. For this ar- 
gument, we work in units where the dimensional- 
ity [energy] = [momentum] = [length] ^"'^= [time] ""'^ and la- 
bel any quantity by the number of powers of energy that 
make it up. For example, energy will have dimension 
and length dimension -1. By noting that the action must 
be dimensionless, the dimensionality of '0 is 1, and A is 
1/2 . According to the standard Wilsonian RG argument 
|28|, only coupling constants with non- negative dimen- 
sion will contribute to the low energy physics. For ex- 
ample, mass terms and the gauge coupling e are relevant 
perturbations. Along with the Chern-Simons term and 
the fermion masses, the terms in the QED3 action ([T]) are 
the only relevant and marginal operators consistent with 
Lorentz and gauge symmetries that can be constructed 
from -0 and A fields. 

IV. 1+1 DIMENSIONAL OPTICAL LATTICE 
QED 

We now discuss the 1-t-l dimensional analog to the 
model discussed in sections |IT] and IIIIl 

The Feynman- gauge Hamiltonian for dimensional 
QED (QED2) is 

H ^Y.^jiivFk + eAi)i^ + ^d^A,d^A, (40) 

k 

where "0 is a two component fermionic spinor, A^, are 
bosonic fields, v — 0, 1, and 71 = Uj,. Following our pre- 
vious arguments one would expect this to be emulated by 
the l-fl dimensional generalizion of our planar Hamilto- 
nian (HH), 
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{x e Lf,y e Lh} 

XX X 

x^Sx + fpR^+i^Rx+l) + H.C. 

b=a.l:k b=0,l\x 

+5l Y '^R/Sx'^R/Sxniy=x±l/2- 
<xy> 



<yx> 

(42) 



The lattice for this model is shown in figure [3 Here, the 
X represent fermion lattice sites (L/), and the y represent 
the boson lattice sites (Lt) associated with every other 
bond between fermion lattice sites, ebk is the boson ki- 
netic term and chemical potential, and it has a minimum 
at K = 71/21. Bosons of species oq are confined to the 
fermion lattice sites and bosons of species Oi occupy the 
y sites. 

Unfortunately, the Hamiltonian in Eq. (PT|) does 
not yield QED2. While this Hamiltonian trivially has 
fermions with a Dirac dispersion, the boson mediated 
fermion hopping term does not lead to a Bose-Fermi cou- 
pling of the form in Eq. Instead, the resulting cou- 
pling, e'0Ai72'i/', has an erroneous extra phase factor of 
e"^/^, and violates Lorentz Invarience. This same phase 



factor was also found in the 2D case, but there it was 
benign: it simply meant that Ax is produced by bosons 
on the y-bonds, and Ay is produced by bosons on the 
x-bonds. In the Id case, however, we cannot make such 
a relabeling (since no 72 matrix appears in the fermion 
kinetic term). We can fix Eq. (PT|) by implementing 
additional phase factors into the mediated hopping ma- 
trix elements. Conceptually, the best way to implement 
these phases is to work in a larger space, allowing the 
fermions to hop onto the y sites, with a total energy cost 
of U relative to the x sites. Using the techniques from 
[21I . m, [3^ , one can engineer a spatially dependent 
phase on the hopping matrix elements, and produce a 
Hamiltonian, 



{x e Lf,y e Lb} 

H = i'R/six±i^s/Rix + ^ly^Ty + 50 51 ^R/Six'^R/Sixnox + 51 XI ^ty^ty'^iy (43) 



Y '^bkUbk + T^bx/y {nbx/y - l) + Y ^ly+2l^V (^Sx'^Sx + 1pRx+li^Rx+l) + H.C. + H, 



2 

6=0,l;fe b=0,l;x <yx> 



mcd 



J^mcd = ^tY e'^^'^ly^R/Slx + H.C. (44) 



<xy> 



Here uJte^^^ is the matrix element for a fermion to hop onto a boson site, with Q ~ 7r/2/, and U is the on-site energy. 
Working to lowest nontrivial order in uit/U, we integrate out the y-sites, obtaining 

H^,, = '§Y [(^k+iV'S.e'^- + 4xi^Rx+,e-^Q-) (1 + ^niy) + i^^/sxi^R/S (l + I"!.)] , (45) 



which gives the desired low energy Hamiltonian. 



V. 3+1 DIMENSIONS 

Full 3+1 dimensional QED can be simulated by com- 
bining the schemes which we have presented for one and 
two dimensions. We will work with a cubic lattice formed 



by layering the 2D square lattice of section HI, with lat- 
tice spacing 1^ between each layer. Each 2D layer is the 
same as the lattice of section HI, containing three bosons 
and one fermion distributed between the various sublat- 
tices. Between the layers we add a fourth boson species 
03, which occupies the space between every other layer 
and can only interact with the fermions through a light- 
assisted transition as in the Id case. We shall then imple- 
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t t 

H,ao} Cf Cf ^ 

Ri Si Ri Si 

FIG. 5: Lattice for the Id theory, with sublattices labeled 
by spin and by the particles that inhabit them. Bold lines 
represent ordinary hopping terms and dashed lines indicate 
hopping facilitated by optical transitions, as described in the 
text. 

ment additional phases to change the sign of the hopping 
amplitudes between planes; we choose the phases such 
that vertical hopping between sites on the R sublattice 
gets no additional phase but z-directional hopping matrix 
elements between sites on the S sublattices are opposite 
in sign. We ignore any vertical hopping beyond near- 
est neighbors. Under those conditions, our free particle 
Hamiltonian is 

H = -2^V^>/7Ak, (46) 

k 

_ f tz COS kzlz t COS kxl + it COS kyl \ 
\. t cos kxl — it cos kyl —tzCoskzlz I 

This Hamiltonian exhibits eight Dirac points, located at 
k = n {±1/21, ±1/21, ±l/2lz}, though many of these are 
degenerate. We can obtain the necessary bispinor rep- 
resentation by simply combining pairs of them; for ex- 
ample, the combination of K = ■k{1/21, 1/21, 1/21^} and 
— K yeilds the following representation for the Lorentz 
group 

The Az component of the gauge field is derived from the 
as bosons in the same manner as in section IIVI Other 
than ensuring that the sign flips between sublattices are 
preserved when adding the gauge interaction, the deriva- 
tion is essentially unchanged. 

VI. SUMMARY AND OUTLOOK 

We have proposed a number of experiments involv- 
ing neutral atoms in optical lattices which will emulate 



the properties of relativistic quantum electrodynamics. 
These models are technically very difficult to realize, re- 
quiring a large number of lasers and atomic species that 
all must be fine-tuned to produce the desired behav- 
ior. However, the necessary technologies (optically in- 
duced phases, interpenetrating spin dependent lattices, 
and tunable Bose-Fermi mixtures) for our model have 
all been separately realized in recent experiments. Fur- 
ther, as a weakly coupled route from a lattice model to 
a Lorentz-invariant dynamical gauge theory, the scheme 
we have proposed here is also of significant theoretical 
interest. 

Finally, it would be interesting to establish whether 
or not this scheme could be extended to simulate more 
complicated gauge theories, such as non-abelian theo- 
ries or pcrturbative gravity. Given the importance of 
these theories in high energy physics, a method for sim- 
ulating them in cold atoms would be a very significant 
development. The difficulty of this process should not 
be underestimated. These theories intrinsically have a 
very large number of degrees of freedom; in 2-1-1 di- 
mensions, an SU(2) theory has nine gauge boson modes 
(three components of spacetime with three SU(2) gen- 
erators associated with each), and an SU(3) theory has 
twenty-four! Further, non-abelian gauge theories have 
boson self-interactions that must be finely tuned to avoid 
breaking gauge invariance. Worse yet, the gauge fields 
exhibit ghostlike degrees of freedom with negative norm. 
Accounting for all of these effects in any scheme would 
be a daunting. Nonetheless, quantum simulation of 
non-abelian dynamical gauge theories is a fascinating 
prospect, and the theory we have described here could 
be an important first step in that direction. 
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